arXiv: 1507.05076vl [astro-ph.GA] 17Jul2015 


Astronomy & Astrophysics manuscript no. ms ©ESO 2015 

July 20, 2015 


Precision age indicators that expioit chemicaiiy pecuiiar stars 

Guy Worthey 

Washington State University, 1245 Webster Hall, Pullman, WA, 99164-2814 
e-mail: gworthey@wsu.edu 

Received 22/05/2015 / Accepted 13/07/2015 


ABSTRACT 

Context. The integrated light of distant star clusters and galaxies can yield information on the stellar formation epochs, chemical 
abundance mixtures, and initial stellar mass functions, and therefore improve our understanding of galaxy evolution. 

Aims. We would like to find a way to improve the determination of galaxy star formation history from integrated light spectroscopy. 
Methods. Several classes of chemically peculiar (CP) stars arise during the course of normal evolution in single stars and noninter¬ 
acting binary stars. An aging stellar population has periods of time in which CP stars contribute to the integrated light, and others in 
which the contributions fade. The HgMn stars, for example, occupy a narrow temperature range of 10500 to 16000 K, which maps to 
a narrow range of ages. Wolf-Rayet stars, He-poor stars, Bp-Ap stars, Am-Fm stars, and C stars all become very common in a normal 
stellar population at various ages between zero and several Gyr, fading in and out in a way that is analogous to features used in stellar 
spectral classification. We examine population fractions and light fractions in order to assess the feasibility of using CP stars as age 
tracers. 

Results. Even though CP stars do not usually dominate in number, there are enough of them so that the CP spectral features are 
detectable in high-quality integrated spectra of young and intermediate age stellar populations. The new technique should be calibrat- 
able and useful. Furthermore, using CP signatures as age dating tools sidesteps reliance on photometry that is susceptible to dust and 
Balmer features that are susceptible to nebular fill-in. 
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1. Introduction 


Tracing stellar population age by integrated light is a mainstay of 
galaxy evolution interpretation. Somewhat ambiguously, yet still 
of great value, blue photometric colors generally indicate galax¬ 
ies with ongoing star formation, while red colors indicate galax¬ 
ies that have not formed significant numbers of stars in the last 
few hundred million years. A technique that often gives higher 
precis ion is using Balmer feature strengths in yo ung popula¬ 
tions dLevesaue & Leithereill201^ or in old ones (iTrager et al.l 
2000l) . While these considerations yield a mean, light-weighted 
age fairly readily, uncertainties remain. First, since all Balmer 
features measure approximately the same thing, namely stellar 
temperature, they are equivalent to one measureable quantity, 
and a mean age is all one can realistically hope to achieve, which 
is rather far from the ultimate goal of uncovering the star forma¬ 
tion history. There are also astrophysi cal degeneracie s such as 
age-metallicity degeneracy at old ages (IWorthey|ll994l) and age- 
initial mass function deg eneracy and age -interstellar medium de¬ 
generacy at young ages (lLeithererll2005l) . 


The motivation for this paper is to highlight a new strategy 
for age dating stellar populations using only high-quality inte¬ 
grated light spectra. Since the method relies on chemically pe¬ 
culiar (CP) stars, a description of the subtypes most likely to be 
useful is in order. 


Wolf-Rayet (WR) stars are mass-losing upper main- 
sequence stars. Because of their extreme winds and subsequent 
mass loss, portions of the WR stars that have undergone nu¬ 
clear processing are exposed and they become CP and exhibit 
emission line features in their spectra. The notable bump in the 
A = 4640 - 4690A region (|C ontilll99l1) is seen in a morpholog¬ 


ically heterogeneous handful of starburstin g galaxies known as 
Wolf-Rayet galaxies (ISchaereret al.lll99(^. This class of galax - 
ies is already being modeled and studied (ILeitherer et al.ll201^ . 

Moving to cooler classes of CP stars, in stars with slow 
rotation and radiative atmospheres a combination of species- 
dependent gravitational settling and radiation pressure can lead 
to heavy species levitating to the top of the photosphere, cau sing 
strong alterations in the observed spectrum (iMichaudI 19701) that 
explain at least some of the large numbers of peculiar spectral 
types among BAF stars (ICowlev et al.lll96^ . We are interested 
in the numerically most common subtypes. They are He-weak 
(CP4) stars that span spectral types B2-B8, HgMn (CP3) stars 
that span spectral types B6-A0, Bp-Ap stars (CP2) that span 
spectral types B6 -F4, and Am -Fm stars (CPI) that span spec¬ 
tral type s A0-F4 (ISmithlll996l) : the designations in parentheses 
are from iPrestonI (Il974t) . 

The diffusion (heavy species levitatio n) timescales are shor t 
compared to the evolutionary timescales (iMichaud et al.l[T97^ . 
and so it is probably a combination of binarism, rotation, 
and magnetic fields that control the levitation CP phenomenon 
(iBailev et al.ll2014t) . The fraction of CP stars as a function of 
age cannot at present be predicted ab initio. However, the stars 
are all main sequence or near main sequence, and are therefore 
amenable to precise empirical scrutiny. 

It is known that t he various subtypes he in well-segregated 
temperature regimes (iNortblll 9931: [Smithll 19961) . which is a cru¬ 
cial aid to their use as extragalactic age indicators. The levitating 
CP stars are common. HgMn stars acco unt for 15% of B8 stars, 
and over 60% of A6 stars are Am or Ap (ISmi thll996l) . For stellar 
populations at the correct ages, their optical-UV spectral features 
should be strongly influenced by chemical peculiarities. 
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Finally, coolest of all, C stars arise late in the lives of 1.5 - 
4 Mq stars on the asymptotic giant branch where a third dredge- 
up conve cts sufficient C to the surface to tip the number ratio 
C/O > 1 (iMarigo et al.ll2008ll . The mass range translates to pop¬ 
ulation ages between 0.5 and 5 Gyr for heavy element abun¬ 
dances typical of the Magellanic Clouds. When C begins to dom¬ 
inate, the spectra of these giants transform from being shaped 
primarily by H 2 O and TiO molecular features to CN, C 2 , CH, 
and similar, while C O remains strong. These stars are very cool 
(iBergeat et al.ll200T!) . but are bolometrically the brightest stars in 
the population. A measure of the ratio of TiO to C 2 features, for 
example, would readily give an estimate of the C star fraction, 
and therefore additional age information. Unlike WR or levitat¬ 
ing CP stars, for which blue or UV spectra are most revealing, 
the telltale spectra for C stars should be in the red or infrared. 

In this Letter, we let the topic of WR stars stand as suffi¬ 
ciently demonstrated, in the sense that it is clear that WR features 
can be seen in the integrated spectra of some galaxies, and that 
seeing such features provides additional leverage on the ages of 
the youngest stars present in the observed galaxy. Section 2 dis¬ 
cusses the feasibility of HgMn stars as age indicators. Section 3 
gauges the observable effect of C stars on intermediate age stellar 
populations, and Sect. 4 briefly discusses the other CP classes. 
Section 5 extends the discussion to the future of this technique, 
outlining the groundwork necessary to make it viable. 

2. HgMn Stars 

According to ISmithI (Il996h . HgMn stars span spectral classes 
AO through B6 and temperatures from 10500 to 16000 K. This 
translates to main-sequence turnoff masses in a range from 3 to 4 
Mq, and that in turn trans lates to a narrow age range of 75 to 150 
Myr at solar metallicity (iBertelli et al.lll99^ . While they exist 
in the stellar population, their num ber fraction a pparently hovers 
near 10% of B7, B8, and B9 stars (ISmithlll996h . 

Their spectral classifications (and names) rest upon a Hg ii 
feature at 43984 and Mn ii features at 43944, 44137, 44206, and 
44282. There are a smatteri ng of HgMn stars in the Coude Feed 
Library dValdes et al.ll2004li . and we investigated example stars 
to get a better idea of the true spectral signature for this CP 
subclass, finding several more strong Mn features. An illustra¬ 
tive pair of CP and normal stars is illust rated in Fig. [H slightly 
smoot hed, along with wavelengths from iKramida & Sansonettil 
(I2013h . 

Also plotted i n Fig. [T| are two average blue galaxy spectra 
from iDobos et al.l ( 1201 2 |k their two bluest categories, SFO and 
SFl. The Mn features are not obvious until the bluest average 
(SFO) is divided by the next (SFl), and then the four strongest 
Mn features seen in the stars also show up as small extra differ¬ 
ential absorption. Most metallic features should be fading as the 
average stellar population becomes hotter, as seen by the indi¬ 
cated Ca and Fe features, and Mn is never strong except in HgMn 
stars, so getting more Mn absorption when moving to hotter stars 
is remarkable. In combination with the fourfold wavelength co¬ 
incidence, this is strong evidence that the HgMn spectral signa¬ 
ture has been detected in aggregate galaxy spectra. It follows that 
the HgMn phenomenon can be used in extragalactic age dating 
once appropriate models are constructed. 

3. C stars 

The stellar evolutionary calculations of iMarigo et alJ (l2008h in¬ 
clude the surface C/O abundance of AGB stars. We included 



Wavelength 

Fig. 1. Portion of the spectrum between and Hy. Av erage galaxy 
spectr a of categories SFO (blue) and SFl (green) from iDobos et all 
(l2012f) and their division scaled about the average by a factor of five 
(black) are shown along with normal B8 sta r HD 152614 ( i Oph ; purple) 
and HgMn star HD 358 ( a And; red! from | Valde s et ^ ll2004l) . Wave¬ 
lengths of Mn transitions dKramida & Sansonettill2oT^ are marked as 
orange vertical lines. The feature of keenest interest is that the four 
strongest Mn features in HD 358 appear as small differential absorp¬ 
tions as one goes from SFl to SFO average galaxies. 

them in the IWorthevI (1 19941) population models, along with syn¬ 
thetic C-rich spectra from lAringer et al.l (l2009h where we im¬ 
posed a correction for self-extinction of Ay = (2800 - 7’efr)/800 
for stars with Tefr < 2800 K in order to compensate for the lack 
of dust i n that particular grid and make it come into agreement 
with the lBergeat et al.l (1200 ll) color temperature relations. 

The low-resolution spectra from these calculations are dis¬ 
played in Fig.|2]for simple stellar populations of solar abundance 
and ages 0.5, 1.0, and 2.0 Gyr, and then displayed again with the 
C stars replaced by M stars of the same atmospheric parameters. 
Redder than R band, the spectral changes are obvious. 

For the purposes of this Letter, we do not adjust the mod¬ 
els for recent suggestions for downward rev isions to the pre¬ 
dicted number of C stars (iGirardi et alJl201^ : a glance at Fig. 
|2] is enough to understand that factors of two or three are not 
sufficient to erase the C star spectra signatures, many of which 
exceed 25% in that figure. Ongoing discussion of the numbers, 
luminosity functions, and frequencies of C stars are healthy steps 
toward using them as age indicators. 

One additional caution about C stars is that they are 
rare, evolved stars susceptible to stochastic counting effects 
that can have effects o n integrated properties in some cases 
(ISantos & Frogeliri997h though entire, large galaxies are prob¬ 
ably immune from stochastic effects. The other, bluer CP stars 
discussed in this Letter are not susceptible, because they are 
main-sequence stars and are therefore very numerous. 

4. Remaining CP Categories 

According to lNorthI (119931) . who studied CP stars that are mem¬ 
bers of open clusters. Si-strong and He-weak stars (also known 
as hot CP2 stars or Bp-Ap stars) are already present at the 5% 
level in clusters of age 10 Myr, though they are not located at 
the main-sequence turnoff (MSTO) at that age, and they hover at 
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Fig. 2. Simple stellar population model spectra at solar abundance and 
three ages, 0.5, 1.0, and 2.0 Gyr, each age labeled on the plot. Car¬ 
bon stars are included (solid lines) or M stars are substituted for C 
stars (dashed lines). While CO absorption stays about the same for both 
cases, there is marked and opposed variation in the H 2 O and TiO fea¬ 
tures of M stars compared to C 2 and CN features due to molecular bal¬ 
ancing. 

10% frequency thereafter until disappearing at around 500 Myr. 
The cooler CP2 stars that show enhanced Sr, Cr, and Eu appear 
at ~60 Myr and disappear at ~1 Gyr, never exceeding 5% fre¬ 
quency. Am stars (CPI) also appear at ~60 Myr at the 5% level 
(again, not at the MSTO), but just as the hot CP2 stars are fad¬ 
ing at 500 Myr, the Am frequency appears to rise to 25% before 
falling off again at or after 1 Gyr. North’s study is limited espe¬ 
cially at the old end by a dearth of clusters. 

These CP categories are shown in context with the others and 
with the rise and fall of the Balmer feature index in Figure 
[3l which is an imprecise amalgamation of data and models. If we 
can properly calibrate all the CP subtypes and then detect their 
presence in integrated light, however, it follows that we can use 
them as precise age dating tools. 

5. Discussion and conciusion 

Considering the comings and goings of hot and cool CP2 stars, 
Am stars, HgMn stars, WR stars, and C stars, we see that some 
event or other happens every few tenths along a log t time¬ 
line (Figure O. This excellent “age resolution” is improved by 
including time-tested photometric colors and Balmer features. 
Given this new wealth of observational signatures, we can realis¬ 
tically consider deriving star formation histories from integrated 
light alone. 

The observability of the CP age indicators is good, though 
high S/N is required for security. To estimate S/N requirements, 
we suppose that a CP star has an increase (or decrease) in 
its various absorptions of equivalent width Wcp- That equiv¬ 
alent width is attenuated by stars in other parts of the HR 
diagram by a factor fpop and by the numerical rarity of the 
CP star category fcp so that the measured equivalent width is 
W - Wcpfpopfcp- The uncertainty depends on S/N as cr{W) - 
V2(A/l-W)/(5 /A), where Ad is the window over which the mea¬ 
surement is made, and the S/N ratio refers to that same spectral 


Fig- 3- Schematic time line showing the age indicators. The index 
dWorthev & Ottavianil 19971) for the hottest star near the main-sequence 
turnoff (MSTO; black broken line) and for the integrated light of the 
simple stellar population of solar abundance (purple dots) refer to the 
right-hand axis. The remainder of the lines refer schematically to the 
left-hand axis. The Wolf-Rayet fraction (blue line) is entirely schematic. 
The lines for He-rich (magenta line), Hg-Mn (slate blue line). Bp (thin 
green line), Ap (dashed orange l ine), and Am (thin black line) stars 
are estimated from data in iNorthI l ll993tl with approximate stellar mass 
converted to a MSTO age and an o verall normal ization estimated from 
Bright Star Catalog classifications dSmithlll9^^ . The line for C stars 
(dash-dotted maroon line) is from the models, an d therefore refers to 
the theoretical predictions of iMarigo et al.l d2Q08h . and is expressed as 
the fraction of C stars compa red with the total number of AGB stars in 
the population. ISmithI d 19961) finds signi ficantly more Am and HgMn 
stars and somewhat fewer Bp stars than iNor t3 dl993h . In this figure, 
the Bp star category includes both Si-strong and He-weak types, which 
transition at a boundary approximately where HgMn stars peak. There 
is probably a sharper boundary between Ap and Bp than is evident in 
this figure, and the old-age Am boundary is uncertain due to scarcity of 
data. 

window (IVollmann & Eversbergll2()()^ . For a single, weak ab¬ 
sorption, if a particular CP type is located at the main-sequence 
turnoff and has 10% frequency {fcp - 0.10) and has a spec¬ 
tral feature that is 95% of continuum and about Ad = 4A wide 
{Wcp ~ 4(0.05) = 0.20 A), and if the turnoff stars contribute half 
the light in the blue {fpop - 0.50), then W =» O.OlA. In that case, 
one would want a continuum S/N ~ 2000 over the wavelength 
span of the weak feature for secure detection. This worst-case 
requirement is greatly lessened in the usual case where many (if 
not hosts of) absorption features contribute to the signal, and also 
many signature absorptions (especially H, He, and Ca) are much 
stronger than 5% deep. 

The S/N requirement also lessens in the UV, where the 
absorptions are stronger and the dilution by non-turnoff stars 
lessens. Working to increase S/N requirements is that CP stars 
often have depressed continuum in the UV, but it still works out 
to be a net positive gain as long as the line absorption increases 
faster than the continuum depression. 

The observability as a function of CP type depends on the 
number and strength of the absorption features in the observed 
spectral window. Roughly speaking, however, relative to Hg-Mn 
stars, all the types we mention will be easier to observe. He-rich 
stars are characterized by a good collection of strong He lines 
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throughout the blue and visual, weak H lines, and all but ab¬ 
sent C lines, the combined equivalent widths of which should 
comfortably overwhelm the Mn lines we point out for HgMn 
stars. Bp and Ap stars are the more extreme cousins of Am stars, 
and show strong Si, Sr, Cr, and Eu, again with a healthy sum 
of equivalent widths over the many features present. Am stars 
have a weak Ca K feature combined with strong metal lines, in¬ 
cluding Fe, and therefore again containing a very healthy sum 
of total absorption. In the Am case, complete disambiguation of 
the Am signature may require measuring Ca K accurately, wor¬ 
risome because it is a single spectral feature. To ease the worry, 
however, we note that this is a strong absorption feature with 
AW ~ 1 to 2A when the Am peculiarity is imposed. 

We did not list some CP categories because of their relative 
rarity. The small number of CP stars lessens their influence on 
integrated light. These include A Boo stars, metal-poor CH stars, 
and stars that are CP because of binary evolution effects such as 
barium stars and Sr4077 stars. These stars should be kept in mind 
for the future, as their spectral signatures could, firstly, slightly 
blur the results from the more common CP varieties, and, sec¬ 
ondly and more positively, some day enable a direct measure¬ 
ment of binary fractions in external galaxies. 

The calibration of the technique is in its infancy and needs 
much more work. Spectral libraries should be built with explicit 
inclusion of sequences of CP stars. If synthetic stellar libraries 
are employed, new grids will have to be computed to accomo¬ 
date the altered abundances, but a more systematic study will be 
needed to specify the abundance patterns that should be synthe¬ 
sized as representative of each class of CP stars. 

Spectral libraries should be extended to the UV, if possible, 
because t he spectral signat ures of chemical peculiarity become 
stronger (lFuhrmannlll989alf5l and the light fraction contributed 
by stars near the main-sequence turnoff increases, though these 
trends are somewhat attenuated by a lowering of the UV contin¬ 
uum in CP stars. Going to the UV is also convenient for applica¬ 
tion to high-redshift galaxies. 

For evolutionary population synthesis, the numbers and tem¬ 
peratures of stars need to be known with confidence as a func¬ 
tion of age and heavy element abundance. In this regard, GAIA 
distances to thousands of nearby CP stars will be very helpful 
regarding temperature spreads and number fractions, but contin¬ 
ued survey work in open clusters, associations, moving groups, 
and well-studied local group galaxies for CP stars is even more 
essential because the ag es and initial metallicities can be known 
dFandstreet et al.ir2007ll . 

In conclusion, this Fetter presents an outline for a technique 
to unravel age structures in the integrated light of stellar popula¬ 
tions using CP stars. 

- The CP technique has high age resolution compared to the 
Balmer feature technique because most CP stars turn on or 
turn off fairly suddenly as a function of population age. 

- The CP technique is insensitive to intervening interstellar 
dust. 

- The CP technique avoids the use of Balmer features, and 
therefore all complications arising from nebular emission. 

- The CP technique relies on high S/N galaxy spectra. 

- The CP technique also relies on a suitably thorough model¬ 
ing effort that includes expanded stellar spectral libraries and 
better empirical constraints on number fractions of CP stars 
as a function of age, chemistry, and evolutionary state. 
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